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Abstract 

We present a new method to explicitly define Abelian functions asso- 
ciated with algebraic curves, for the purpose of finding bases for the rele- 
vant vector spaces of such functions. We demonstrate the procedure with 
the functions associated with a trigonal curve of genus four. The main 
motivation for the construction of such bases is that it allows system- 
atic methods for the derivation of the addition formulae and differential 
equations satisfied by the functions. We present a new 3-term 2-variable 
addition formulae and a complete set of differential equations to generalise 
the classic Weierstrass identities for the case of the trigonal curve of genus 
four. 
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1 Introduction 



In this paper we present a new addition formula and sets of differential equations 
satisfied by the Abelian functions, with poles along the standard theta divisor, 
associated with a trigonal curve of genus four. Perhaps of greater importance, 
is the discussion of how similar results may be derived for any suitable curve 
through the derivation of bases for the corresponding vector spaces. The prob- 
lem of choosing a basis of linearly independent functions from a wider set can 
be solved using the a- function expansion, (see for example [12] and [18]). The 
problem of identifying the sets of suitable functions from which to choose is the 
main focus of this paper. As part of the discussion, we present a number of new 
explicit classes of such functions which may be used for any suitable curve, (i.e. 
without restriction on the degree of the curve variables or the genus). 

We apply these new results to the trigonal curve of genus four, allowing us 
to derive a new addition formula in Theorem [8] The existence of the addition 
formula was shown in [5] but the formula could only now be derived after the 
explicit derivation of a basis for the vector space of Abelian functions having 
poles of order at most three along the standard theta divisor. A similar addition 
formula has recently been derived for the trigonal curve of genus three in |15j . 
This also required the derivation of a new basis which used the general ideas 
discussed here. Such formulae were also the topic of [14] which presented results 
for the Weierstrass p and cr-functions along with generalisations to genus two. 

We also derive complete sets of differential equations for the cyclic trigonal 
curve of genus four. This can be seen as a continuation of the work started in 
[5 , however the complete sets presented here required both the knowledge of 
the bases and new efficient computational techniques first introduced in |18j . 
Some of the relations are very lengthly and so not included in the paper itself. 
However, they are stored in a variety of formats online at [16) . 

Many of the results we present can be viewed as generalisations of classic 
results for elliptic functions. We will conclude the introduction below by remind- 
ing the reader of the relevant results from the theory of Weierstrass functions. 
Then in Section [2] we revise the definitions of the Abelian functions, discussing 
the general properties they satisfy. We proceed in Section [3] to consider the 
problem of determining bases for the vector spaces of such functions, presenting 
an explicit construction for the trigonal curve of genus four. In Section [4] we 
derive sets of differential equations satisfied by the functions and in Section [5] 
we present the new addition formula. 

The classical results for Weierstrass elliptic functions forms a template for 
our theory. Let p(u) be the Weierstrass p-function, which as an elliptic function 
has two complex periods u\ , u>2- 

p(u + uji) = p(u + uj 2 ) = p(u), for all u 6 C. (1) 

The p-function has the simplest possible pole structure for an elliptic function 
and satisfies a number of interesting properties. For example, it can be used to 
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parametrise an elliptic curve, 



y 2 = 4x 3 -g 2 x-g 3 , (2) 

where g 2 and 173 are constants. It also satisfies the following well-known differ- 
ential equations, 

(p'{u)) 2 = 4p(«) 8 - g 2 p(u) - S3, (3) 
p» = &p(u) 2 -\g 2 . (4) 

Weierstrass introduced an auxiliary function, <r(u), in his theory which satisfies 
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p( u ) = -^2 lo § ( 5 ) 

The cr-function plays a crucial role in both the generalisation and in applications 
of the theory. It satisfies the following two term addition formula, 

a(u + v)a(u — v) , . . . 

In this document we present generalisations of equations (JTJ) — © for higher 
genus functions. 

Elliptic functions have been the subject of much study since their discov- 
ery and have been extensively used to enumerate solutions of non-linear wave 
equations. Recent times have seen a revival of interest in the theory of Abelian 
functions, which have multiple independent periods, and so generalise the ellip- 
tic functions. The periodicity property is usually defined in association with an 
underlying algebraic curve. These functions have been shown to solve some of 
the differential equations arising in mathematical physics and have been used 
in a variety of applications. 



2 Constructing Abelian functions 

In this paper we study Abelian functions associated with the following set of 
algebraic curves. 

Definition 1. For two coprime integers (n, s) with s > n we define an (n,s)- 
curve as a non-singular algebraic curve defined by f(x,y) — 0, where 

f(x,y) = y n +pi{x)y n - 1 + p 2 {x)y n - 2 + ■ ■ ■ +p n -i(x)y - p n {x). (7) 

Here x, y are complex variables and Pj{x) are polynomials in x of degree (at 
most) \Js/n\. We define a simple subclass of the curves by setting Pj(x) = 
for < j < n — 1. These curves are defined by 

f(x, y)=y n - (x s + Xs^x 3 - 1 + ■ ■ ■ + \ x x + A ) (8) 
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and are called cyclic (n, s) -curves. We follow previous publications and de- 
note the curve constants by Xj for the cyclic curves and by Uj for the general 
(n, s)-curves. Note that in the literature the word "cyclic" is sometimes replaced 
by "strictly" or "purely (n-gonal)" . 

Unless we specify otherwise, we will be working with the general (n, s)- 
curves. In some cases we restrict to the cyclic case for either theoretical or 
computational reasons, as noted at the time. 

We denote the surface defined by an (n, s)-curve by C. The genus of C is 
given by g = \(n — l)(s — 1) and the associated Abelian functions, to be defined 
shortly, will be multivariate with g variables, u = (ui, . . . , u g ). As an example, 
the elliptic curve in equation J2]) is a (2,3)-curve and the associated Weierstrass 
a and p- functions depend upon a single complex variable u. 

The (n, s)-curves with n = 2 are generally defined to be hyperelliptic curves, 
(when s = 3 they reduce to elliptic curves). Klein developed an approach to 
generalise the Weierstrass p-function to Abelian functions associated with hy- 
perelliptic curves, as described in Baker's classic texts [2] and [4]. This approach 
has motivated the general definitions in [7J and [T3] of what we now call Kleinian 
p-functions. It is the properties of these and the generalised cr-function that are 
our objects of study. 

In the last few years a good deal of progress has been made on the theory 
of Abelian functions associated to those (n, s)-curve with n — 3, which we 
label trigonal curves. The functions associated to a trigonal curve were first 
studied in detail in [5], where the authors of [7J furthered their methods for 
the hyperelliptic case to obtain realisations of the Jacobian variety and some 
key differential equations between the functions. The simplest trigonal curve 
is the (3,4)-curve which is of genus three. This was first studied in detail in 
[12] with new results recently presented in [15]. The other canonical class of 
trigonal curve is the (3,5)-curve, which has genus four and will be the example 
considered in this paper. This class of curves is defined by 

y 3 + (p 3 x + p 5 )y 2 + (n lX 3 + piiX 2 + [i 7 x + pi W )y 

= X 5 + fl 3 X 4 + pL 6 X 3 + fJ, 9 X 2 + pi 12 X + ^15, (9) 

with the cyclic restriction given by 

y 3 = x 5 + X 4 x 4 + X 3 x 3 + X 2 x 2 + Aix + A . (10) 

The Abelian functions associated with such curves were first studied in [B] for use 
in an application to the Benney equations, shortly followed by the more detailed 
study of the curves given in [S] . The authors of these papers provided an explicit 
construction of the differentials on the curve, solved the Jacobi inversion problem 
and obtained a number of differential equations between the p-functions. Some 
of the properties of higher genus trigonal curves have recently been explored in 
[17] . The class of (n, s)-curves with n — 4 are defined as tetragonal curves and 
were considered for the first time in [18] and [19]. The lowest genus tetragonal 
curve is of genus six, leading to some computational restrictions. 
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We define a set of weights for the theory, denoted by wt and often referred 
to as the Sato weights. We start by setting 

wt(x) = —n, wt(y) = — s (11) 

and then choosing the weights of the curve parameters to be such that the 
curve equation is homogeneous. We see that for cyclic curves this imposes 
wt(Aj) = —n(s — j) while for the non-cyclic curves we usually label the fJLj with 
their weight. For example, in equation ^ the [ij all have weight —j and the 
equation has weight —15 overall. The weights of all other objects in the theory 
may be derived uniquely and will be commented on throughout the paper. Note 
that all the objects in this paper have a definite Sato weight and all the equalities 
are homogeneous with respect to this weight. A more detailed discussion of the 
weight properties may be found elsewhere, for example in [18] . 

We now discuss how to construct the functions associated to the curve. We 
start by choosing a basis for the space of differential forms of the first kind. 
These are the differential 1-forms which are holomorphic on C. There is a 
standard procedure to construct this basis for an (n, s)-curve, (see for example 
[2Tj). Any (3,5)-curve has a basis given by 

/ dx xdx ydx x 2 dx \ 
^fy(x,y)' fy(x,y)' f y (x,y)' f y {x,y)J' 

In general the basis is given by hidxj f y , i = 1 . . . g where the hi are monomials 
in (x,y) whose structure may be predicted by the Weierstrass gap sequence. 
We note here that some authors use a different normalisation when working 
with (n, s)-curves. For example, in the hyperelliptic cases it is common for 
the coefficient of x s to be set to 4 instead of 1 in the curve definition. This 
corresponds to the removal of a factor of ^ from the basis of differentials, which 
would otherwise each have denominator f y = 2y. It is simple to move between 
such different normalisations. 

We next choose a symplectic basis in Hi(C, Z) of cycles (closed paths) upon 
the compact Riemann surface defined by C. We denote these by {c*i, . . .a g , 
/3i, . . . (3 g }, with g = 4 in the (3,5)-case. We ensure the cycles have intersection 
numbers 

en ■ ctj = 0, Pi - (3j =0, at ■ j3j = 5 l} = 

The choice of these cycles is not unique, but the functions we work with will be 
independent of the choice. 

We introduce dr as a basis of differentials of the second kind. These are 
meromorphic differentials on C which have their only pole at oo. This basis is 
usually derived alongside a fundamental differential of the second kind which 
plays an important role in the theory of p- functions. We do not explicitly use 
either of these objects in this paper and so refer the reader to [21] for the general 
theory and [6] which contains all the details for the trigonal curve of genus four. 



f 1 if i = i, 
\ if i ^ j. 
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We can now define the standard period matrices associated to the curve as 

w '=(^L =1 , w "=(^^) M=1 ,...,, 



"' U : > 'H,, , s ^0M r V M=1 , 

We define the period lattice A formed from lj',lu" by 

A = {uj'm + ui"n, m,neZ 9 }. 

Note the comparison with equation ([T]) and that the period matrices play the 
role of the scalar periods in the elliptic case. All the functions we treat in this 
paper are defined upon C 9 with coordinates usually expressed as 

u= (ui,...,u g ). (13) 

Note that any point ueC 9 can be expressed as 

where the Pi are points upon C. The period lattice A is a lattice in the space 
C 9 . Then the Jacobian variety of C is presented by C 9 /A, and is denoted by J. 
We define k as the modulo A map, 

k : C 9 — > J. (14) 

For k = 1, 2, ... define 21, the Abel map from the fc-th symmetric product 
Sym fc (C) of C to J by 

21 : Sym fe (C) 

(P 1 ,...,P k ) i ^ 1/ du+-+ I du\ (mod A), (15) 

where the are again points upon C . Denote the image of the fc-th Abel map 
by and define the k-th standard theta subset (often referred to as the fc-th 
strata) by 

QlM = W [k] U [-l]iy [fcl , 

where [—1] means that 

•••,%) = {-ui,...,-u g ). 

When k — 1 the Abel map gives an embedding of the curve C upon which we 
define £ = x~ ~ as the local parameter at the origin, 2ti(oo). We can then express 
y and the basis of differentials using £ and integrate to give series expansions for 
u. We can check the weights of u from these expansions and see that they are 
prescribed by the Weierstrass gap sequence. In particular, for the (3,5)-curve 
the variables u — (ui, . . . , U4) have weights (7, 4, 2, 1). 

We now consider functions that are periodic with respect to the lattice A. 




G 



Definition 2. Let dJl(u) be a meromorphic function of u G C 9 . TTien 971 is 
a standard Abelian function associated with C if it has poles only along 
K -i(Q[9-i]) an d satisfies, for all £ G A, 

m(u + £) = £Dl(u). (16) 

We define generalisations of the Weierstrass functions which satisfy equation 
(|16[) . defined using the quasi-periodic function defined below. First, let S = 
oj'S' + ui"6" be the Riemann constant with base point oo. Then [S] is the theta 
characteristic representing the Riemann constant for the curve C with respect 
to the base point oo and generators {ctj, f3j} of i?i(C, Z). (See for example [7J 
pp23-24.) 

Definition 3. The Kleinian cr-function associated to a general (n, s)-curve 
is defined using a multivariate 9-function with characteristic d as 

tr(u) = cexp(iM7 7 Vr 1 M T ) ■ 0[6}((u')- 1 u T | (w')-V). 

2iri{ 

i(m + 6 , ) T (u')- 1 u , \m + *') + {m + <5 / ) T ((w / )- 1 w T + «")} ■ 

XTie constant c is dependent upon the curve parameters and the basis of cycles 
and is fixed later, following Lemma [2j 

We now summarise the key properties of the cr-function. See [7] or [21] for 
the construction of the cr-function to satisfy these properties. For any point 
u G C 9 we denote by u' and it" the vectors in R 9 such that 

u = J -a' + J'u". (17) 
Therefore a point I G A is written as 

£ = w'£' + w"£" G A, G Z 9 . (18) 

For u,v <eC 9 and £ G A, define L(u,v) and x(^) as follows: 
£(«,«) = u T (r?V + ?7'V'), 

X (£) = e Xp [27ri{(£') T ^-(O T ^ + |(0 T ^"}]- 
Lemma 1. Consider the cr-function associated to an (n, s)-curve. 

• Lt is an entire function on C 9 . 

• Lt has zeros of order one along the set k^ 1 (Q^ 9 ^ 1 ^). Further, we have 
cr(tt) outside the set. 

• For all u G C 9 ,£ G A the function has the quasi-periodicity property: 

a(u + £) = x(£)exp 



cexp 



— 1 T 



exp 



2'* 



cr(w). 



(19) 
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• It has definite parity given by 



a(-u) = (-l)n(« 2 -i)(s 2 -iV( M ). (20) 

Hence the function is even in the (3,5)-case. 

Proof. The function is clearly entire from the definition, while the zeros and 
the quasi-periodicity are classical results, (see [2]), that are fundamental to 
the definition of the function. They both follow from the properties of the 
multivariate (^-function. The parity property is given by Proposition 4(iv) in 

□ 

We next define p- functions using an analogy of equation ([5]). Since there is 
more than one variable we need to be clear which we differentiate with respect 
to. We actually define multiple p-functions and introduce the following index 
notation. 

Definition 4. Define m-index Kleinian p-functions for m > 2 by 

where i\ < ■ ■ ■ < i m G {1, . . . ,g}. 

The m-index p-functions are meromorphic with poles of order m when 
er(it) = 0. We can check that they satisfy equation (|T51) and hence they are 
Abclian. The m-index p-functions have definite parity with respect to the 
change of variables u — > [— l]u. This is independent of the underlying curve, 
with the functions odd if m is odd and even if m is even. Note that the ordering 
of the indices is irrelevant and so for simplicity we always order in ascending 
value. 

When the (n, s)-curve is chosen to be the classic elliptic curve then the 
Kleinian cr-function coincides with the classic cr-function and the sole 2-index 
p-function coincides with the Weierstrass p-function. The only difference is the 
notation with 

p n {u) = p(u), p nx (u) = p'(u), pimO) = p"{u). 

Clearly, as the genus of the curve increases so do the number of associated 
p-functions. In general, the number of m-index p-functions associated with a 
genus g curve is 

(g + m — 1)1 
m\(g-l)\ ' 

For the (3,5)-curve we study we have g = 4 giving ten 2-index p-functions, 
twenty 3-index p-functions etc. By considering Definition [4] we see that the 
weight of the p-functions is the negative of the sum of the weights of the vari- 
ables indicated by the indices. We note that curves of the same genus will, 
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notationally, have the same p-functions. However they exhibit different be- 
haviour as indicated by the different weights. When considering curves of genus 
four, we have the trigonal curve which we study and also a hyperelliptic curve, 
the (2,9)-case. The weights of the associated 2-index p-functions are compared 
below. 



Curve 


Pn 


Pl2 


Pl3 


Pu 


P22 


P23 


P24 


P33 


P34 


P44 


(2,9) 


-14 


-12 


-10 


-8 


-10 


-8 


-6 


-6 


-4 


-2 


(3,5) 


-14 


-11 


-9 


-8 


-8 


-6 


-5 


-4 


-3 


-2 



Wc now introduce a final result detailing how the functions can be expressed 
using series expansions. 

Lemma 2. The Taylor series expansion ofo~(u) about the origin may be written 
as 

oo 

a(u)=K-SW„ lS {u) + J2 C k(u,X). (21) 

k=0 

Here K is a constant, SW n , s the Schur-Weierstrass polynomial generated by 
(n,s) and each Ck a finite, polynomial composed of products of monomials in u 
of weight k multiplied by monomials in the curve parameters of weight —(wt(o~) — 
k). 

Proof. We refer the reader to [21] for a proof of the relationship between the 
er-function and the Schur-Weierstrass polynomials and note that this was first 
discussed in [8] . We see that the remainder of the expansion must depend on 
the curve parameters and split it up into the different Ck using the weight 
properties. We can see that each Ck is finite since the number of possible terms 
with the prescribed weight properties if hnitc. 

□ 

Large expansions of this type were first introduced in [Bj in relation to the 
study of reductions of the Benney equations. Since then they have been an 
integral tool in the investigation of Abelian functions. Recently computational 
techniques based on the weight properties have been used to derive much larger 
expansions and we refer the reader to [18] and [17] for details of how to construct 
and use these expansions. We note that such expansions are possible for the 
general curves, but that the calculations involved are far simpler for the cyclic 
cases. 

The Schur-Weierstrass polynomials are Schur polynomials generated by a 
partition associated with (n,s), (see [8] for full details on the construction). 
The Schur-Weierstrass polynomial generated by (3,5) is 

SW 3t5 = 4jgW| - \u\u\ - jU^ + ul+ U Z u\u 2 - U4"l- (22) 

In Definition [3] we fix c to be the value that makes the constant K = 1 in the 
above lemma. Some other authors working in this area may use a different 
constant and in general these choices are not equivalent. However, the constant 
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can be seen to cancel in the definition of the p-functions, leaving results between 
the functions independent of c. Note that this choice of c ensures that the 
Kleinian cr-function matches the Weierstrass cr-function when the (n, s)-curve is 
chosen to be the classic elliptic curve. 

The connection with the Schur- Weierstrass polynomials also allows us to 
determine the weight of the cr-function as (l/24)(n 2 — l)(s 2 — 1). In the (3,5)- 
case this gives cr(u) weight 8 . The cr-function expansion has been calculated 
up to and including the polynomial C50 in the (3,5)-case. These expansions are 
stored online at [IB] . 

3 Bases for the vector spaces of Abelian func- 
tions 

In this section we discuss bases for the vector spaces of Abelian functions. We 
first identify the vector spaces considered, then define some new generic Abelian 
functions before presenting some explicit results for the (3,5)-curve. We note 
that everything in this section is valid for general (n, s)-curves rather than the 
cyclic restriction used in later sections. 

3.1 Classifying Abelian functions by pole structure 

We can classify the Abelian functions according to their pole structure. It is 
usual to denote by r( J, 0(m0' fe l)) the vector space of Abelian functions defined 
upon J which have poles of order at most m, occurring only on the fcth strata, 
0[ fe l. The case where k = g — 1 is of interest since all the Abelian functions we 
deal with have poles occurring here, on the 0-divisor. Hence, for brevity, we just 
refer to such spaces as T(m) throughout the rest of this paper. A key problem 
in the theory of Abelian functions is the generation of bases for these vector 
spaces. Note that the dimension of the space T(m) is m 9 by the Rlemann-Roch 
theorem for Abelian varieties, (see for example )20j). 

Recall that the m-index p-functions all have poles of order m and so are 
natural candidates for the construction of bases for these spaces. In fact, for 
a genus one curve they are sufficient since the Weierstrass p-function and its 
derivatives can span all such spaces: 



In general, the number of p-functions with poles of order at most m is 



r(o) 
r(2) 
r(3) 
r(4) 



r(i) = ci, 

Cl © Cp, 
CI © Cp © Cp', 
Cl © Cp © Cp' © Cp' 



ft 



m 



(g + N-l)\ = (g + m)\ 
Nl(g-1)\ glm\ 



E 



(.9 + 1), 



N=2 
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which, as g increases, will not grow as fast as the dimension of the space. We 
hence need to identify a wider class of Abelian functions than the p-functions 
in order to construct such bases. Note that the p-functions are not sufficient 
even for g = 2. In this case 



but an additional function must be included alongside the {pijk} in the basis for 
r(3). This space has dimension 9 while there are only four 3-index p-functions 
to add. The final entry is filled by the function p\xp%% — p\ 2 - Although each of 
these terms has poles of order four, together they cancel to leave poles of order 
three. The bases for subsequent T(m) in the g — 2 case can be constructed 
from the entries in the basis for T(m — 1) together with the derivatives of those 
entries. This may be verified by noting that the dimension of the spaces and 
the number of distinct function obtained through differentiation increase by the 
same amount each time. It can also be checked that these function all have 
distinct weights and so must be linearly independent. 

These cases of genus 1 and 2 curves, where the entries for the bases may 
be determined in general, are special. They fall into the class where the theta 
divisor is non-singular and the 2?-module structure of such cases is discussed in 
[10] . However, subsequent (n, s)-curves are not covered by this case and so new 
methods must be used to derive bases. One of the authors of [TU] has conjectured 
the general hyperelliptic 2?- module structure in [23], and explicit results have 
been derived for the genus 3 and 4 hyperelliptic curves in [15] . However, this 
paper discusses approaches relevant for all (n, s)-curves, including those that 
are non-hyperelliptic such as our example curve, the (3,5)-case. 

Recall that an entire Abelian function must be a constant and that there 
is no Abelian function with a single pole of order one. Hence those Abelian 
functions with poles of order two are the simplest and so are often referred to as 
fundamental Abelian functions. The basis problem has been solved in general 
for such functions, through the inclusion of the following extra class of Abelian 
functions. 

Definition 5. We define an operator Tj, now known as Hirota's bilinear oper- 
ator although it was used much earlier by Baker in J^j. 



Then an alternative, equivalent definition of the 2-index p-functions is given by 



r(o) 
r(2) 



r(i) = ci, 

Cl © Cpn © Cpi2 © Cp 22 , 



Pij{u) = - TiTj(j{u)(j{v) i < j e {l,...,g}. 

Z(7[U) v—u 



(24) 



We extend this to define the m-index Q-functions, for m even, by 



Qi u i2,...,i m (u) = - Ti Ti ...T lm a(u)a(v) 

Lo~\uY v 




(25) 
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where ii < ■ ■ ■ < i m £ {1, . . . , g}. 

The m-index Q-functions are Abelian functions with poles of order two oc- 
curring when o~(u) = 0. Note that if you were to apply the definition with to 
odd then the functions would be identically zero. They are a generalisation of a 
function hrst used by Baker, referred to as the Q-function. The generic 4- index 
Q-functions were introduced when research first started on the trigonal curves 
(and in the literature are also sometimes just referred to as Q -functions) . The 
definition above was developed in [18] as increasing classes are required to deal 
with cases of higher genus, (see also |T7]). In this paper we only need to use 
4-index Q-functions, which in |12) were shown to satisfy, 



Similar expressions have been found for the higher index Q-functions, (see [18]). 
Note that in some specific cases, alternatives to the Q-functions are available. 
For example, in the (2,7)-case the basis for T(2) requires one additional entry 
after including the p-functions. The general approach would lead to the addition 
of a 4-index Q-function, Q1333, but the function, 



may be used instead. This A-function, originally introduced by Baker in [5], has 
poles of order at most three in general. (This may be checked using Definition 
But in the case of the (2,7)-curve, these can be shown to cancel to leave 
poles of order two. 

The use of A can be advantageous since it allows the theory to be completely 
realised in terms of 2 and 3-index p-functions, with all higher index p-functions 
given recursively in terms of these. Similar functions have been found for the 
genus four hyperelliptic curve, (see [15]), but they appear to be a feature unique 
to the hyperelliptic cases. We note that the A-function also arises naturally 
in the representation theoretic approach to hyperelliptic Abelian functions de- 
scribed in [J. In this paper, we concentrate on deriving classes of functions 
whose maximal pole order is known for any underlying curve, (such as the Q- 
functions), as opposed to searching for functions such as A which are only of 
use to problems on specific curves. The approach presented here may hence be 
used in any applications which require the derivations of such bases. 

When constructing a basis for F(m) we can start by including the entries 
in the basis for T(m — 1). We then know that the remaining entries have poles 
of order to. A natural place to look for these functions is in the derivatives of 
the functions in T(m — 1). Note that while the derivatives of (to — l)-index 
p-functions are m-index p-functions, the same is not true for the Q-functions. 
For brevity we adopt the notation 



Qijke = pijke — Zpijpke — 2pikpje — ^Pupjk- 



(26) 



A = pnp33 - P12P23 - Pl3 + P13P22- 



(27) 



d m Qijkl (u) 



d 



Qijki(u) 



du 



■m 
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and similarly for other functions. 

As discussed in [22], the derivatives of existing basis functions will not be 
sufficient to find successive bases in the case where the theta divisor has singular 
points. In Section 8 of [T2] the authors introduced a new class of functions, p^\ 
defined to be the minor of the matrix [py]3 X 3- The authors used these 

functions to overcome the problem for the three pole basis in the (3,4)-case and 
they appeared again for use in the (2,7)-case in [15]. We can check that these 
functions have poles of order at most three by substituting the p-functions for 
their definition. This motivates a more general definition for similar functions, 
valid for curves of any genus. 

Definition 6. Consider the matrix [pijjgxg- Denote pl J J' fe ''l to be the determi- 
nant of the submatrix formed by rows i and j and columns k and I. (Note that 
we need i ^ j and k ^ I). So 



p[i,j,k,l] 



Pik Pii 
Pjk Pjl 
Pikpjl - Pilpjk- 



We refer to these as cross product p-functions. 

So the genus three p^'-functions used in [T2] are the cross product functions 
in the genus three case. 

Note that while each of the terms has poles of order four, together they 
cancel to leave poles of order three. To see this substitute using Definition 2] so 
that 

d d r i OiO-j - croij 
P ^-c^ l0g ^=— ^ 

and so 

[i,j,k,l] _ ~°~ikV3°~l ~ 0~]lO~iO~k + O'O'ik^jl + O'ijO'kO'l + O ' jlO~ k ~ 0~Q~ijO~kl 
P 5 I 

0~ 

from which it is clear that poles will be at most order three. Although this class 
of pl 1 -^''] functions completed the basis for functions with poles of order at 
most three in the (3,4)-case, it was not sufficient to complete the corresponding 
basis in the (2,7)-case. These two genus three curves were recently discussed in 
detail in [TS]. The difference between the bases resulted from the fact that two 
of the cross product functions were linearly dependent in the hyperelliptic case 
while they were independent in the trigonal case. This is made possible through 
the different weight structures in the two cases. In [TS] the final entry in the 
(2,7) basis for T(3) was filled by the function 

T = P222 + 2 P22 - P22P2222, 

which although constructed from terms with poles of order six, has poles of order 
three overall. This particular function belonged to a wider class of T-functions 
derived through an attempt to match, in general, the poles of a quadratic term 
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in the 3-index p-functions with a polynomial in fundamental Abelian functions. 
These T-functions have poles of order three for any (n, s)-curve. The discovery 
of this class prompted the systematic study of how the p-functions may be com- 
bined to produce classes of functions with prescribed orders of poles, discussed 
in Section 

The classes of functions discussed in the next section are just the candidates 
for the basis entries. It is of course still necessary to check which individual 
functions are linearly independent. This may be achieved for particular cases 
using the series expansion for the cr-function. The calculations involved can be 
lengthly and grow in CPU time and memory requirement with both the genus 
and the number of poles. Significant computational simplifications can be made 
by writing procedures that take advantage of the weight structure present in the 
theory. For example, when expanding the product of series it is only necessary 
to multiply those terms which will give the correct final weight, as the other 
terms must all cancel. Many of the results in this paper required a new efficient 
code to evaluate arbitrary products of series at appropriate weights. The basis 
in Theorem [3J the addition formula in Theorem [5] and the differential equations 
in Section 2] all made use of this and the code is available online at [TB] . 

We can further simplify the basis calculations by noting that there is a finite 
weight range that basis entries can take. Entries in T(m) can have weight no 
lower than — mwt(er), (see Lemma 3.4 in |15j). This stops us from testing those 
functions at a lower weight thus drastically reducing the amount of computation 
required to find such bases. These lower weight functions must be expressible 
as a linear combination of the basis functions in which every term depends on 
the curve parameters. We note that while these bases are not unique, (they can 
formulated using different functions), their weight structures are. We conjecture 
that each basis must have an element with the minimal weight and note that 
the structure of these bases are still under investigation. 

3.2 Classes of functions for generic curves 

The new classes of functions presented in this section were derived in turn 
by considering a prescribed order of poles and a set number of indices on the 
functions. A general polynomial was constructed from p-functions in which 
each term has the set number of indices. The arbitrary coefficients in this 
polynomial were then derived by substituting for Definition U and ensuring the 
poles of higher order are canceled. If all coefficients need to be set to zero then 
no new class is discovered, however, there is usually a non-zero combination 
that gives the required pole order. 

Since the basis problem for T(2) may be solved in general by the Q-functions, 
we start by considering functions with poles of order three. When restricted to 
two indices the only polynomial that may be constructed is a linear combination 
of 2-index p-functions which will have poles of order at most two. Similarly, 
when restricted to three indices we can only form a linear combination of 3-index 
p-functions. The first non trivial case occurs when we allow four indices. We 
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may construct polynomials from 4-index p-functions and the products of two 
2-index p-functions. There are two non-zero combinations for which the poles 
of order four cancel. We identify these by considering an arbitrary sum of pijki 
and the different products of 2-index p-functions and imposing restrictions on 
the coefficients so that the poles of order four cancel. The two combinations 
turn out to be the cross-product functions from Definition [5] and the 2-index 
Q-functions from Definition [5l (in which the poles of order three also cancel). 

When we restrict to five indices then we may construct polynomials from 
5-index p-functions and the products of a 2 and 3-index p-function. There 
are two non-zero combinations for which the poles of order 5 and 4 cancel. 
Again, these were identified by considering an arbitrary sum of pijkik an d the 
different products of 2 and 3-index p-functions and imposing restrictions on the 
coefficients so that the poles of order four and five cancel. The first combination 
is given by 

Bijklm = PijPklm + ^(pjkpilm + PjlPikm + PjmPikl 

— 2pklPijm — ZpkmPijl ~ ^PlmPijk) , (28) 

while the second can be shown to be a combination of £>-functions and deriva- 
tives of 4-index Q-functions. The label B is used as they are bilinear in the 2 
and 3-index p-functions which gives them a special role in the construction of 
differential equations between the p-functions, discussed further in Section FOl 
Note that the order of the indices in the p and Q-functions is irrelevant and so 
we always write them in ascending order. However, this is not the case for the 
£>-functions and so care must be taken when choosing indices. (It is possible to 
replace the ^-function by four separate functions which only need be evaluated 
in ascending indices.) 

We next allow six indices meaning we can build polynomials from the set of 
functions 

{pijklmm PijkPlmn: Pij Pklmn i PijPklPmn\ • 

There are three independent non-zero combinations in which the poles of order 
6,5 and 4 cancel. First there are the 6-index Q-functions as given by Definition 
[5l (in which the poles of order 3 cancel also) . Secondly there is the class 

Tijklmn = Pijkpl ran PkmPln + ^PijPknPlm 

~\~^PikPjlpTnn ~t~ ^PikPjmPln ~t~ ^PikPjnPlm ~t~ ^PilPjkPmn 
+4:pimPjkPln + ^PinPjkPlm ~ ^PilPjmPkn ~ SpilPjnPkm 
SpimPjlPkn - SpimPjnPkl ~ ^PinPjlPkm ~ ^PinPjmPkl 

^Pij Pklmn ^PikPjlmn ^PjkPilmn ^PlmPijkn ^PlnPijkm 

-2p 

mnpijkl ~t~ pilpjkmn ~t~ pimpjkln pinpjklm ~t" pjlpikmn 
+ PjmPikln + PjnPiklm + Pklpijmn + PkmPijln + PknPijlm) 

which are the same as the ^-functions in [15] with the Q-functions replaced by 



15 



their expression in p-functions. Finally, there is a class given by 

^ijklmn Pijklmn ^{pijkPlmn ~t~ PijlPkmn ~t~ PijmPkln ~t~ PijnPklm 

+ PiklPj mn PikmPjln T" piknPjlm T" pilmpjkn T" pilnPjkm T" pimnPjklj- 

This process may be continued indefinitely by increasing the number of indices 
allowed. The process of determining coefficients so that higher order pole cancel 
is not very computationally intensive as it just involves deriving and solving sets 
of linear equations. However, the process of checking whether new functions 
can be formed from derivatives of existing ones can become cumbersome. The 
functions defined here, along with some others that were not used in this paper, 
can all be found online at [16]. The only other class that is used explicitly in 
this paper is given below. 

'Pijklmno — PijPklPmno PijPkmPlno ^PijPknPlmo Pij PhnPkno 

13 1 

+ 2pijPlnPkmo + 2pijPnoPklm + jPijP 

mn Pklo Pik Pjl Pmno 

— pikPjmplno + \pikPjnplmo ~ Pikplmpjno + \pikplnPjmo 

1 3 

+ 2pikPmnPjlo + ^PikPnoPjlm ~ PilPjkPmno ~ PilPjmPkno 

+ 2pHPjnPkmo ~ PilPkmPjno + 2 Pil Pkn Pjmo + 2PHP 

rnnpjko 

3 1 
+ 2pilPnoPjkm ~ PimPjkPlno ~ PimPjlPkno + ^PimPjnPklo 

113 

— pimPklPjno + 2pi™-PknPjlo + ^PimPlnPjko + ^PimPnoPjkl 

3 

— pinPjkPlmo ~ PinPjlPkmo ~ PinPjmPklo ~ ^PinPjoPklm 

3 

— pinPklPjmo — PinPkmPjlo ~ 2~PinPkoPjlm ~ PinPlmPjko 
3 3 

2pinPloPjkm "^PinPmoPjkl PjkPloPimn PjkPmoPiln 
+ Pjkpnopilm + Pjlpkopimn H~ Pjlpmopikn H~ Pjlpnopikm 
J rpjmpkopiln + Pjmpnopikl + Pjmplopikn ~ \pjnpkopilm 

"^PjnPloPikin ^PjnPmoPikl PjoPklPimn PjoPkmPiln 

— \pjopknpilm + Pjoplmpikn ~ \pjoplnpikm ~ \pjopmnpikl 

1 1 

+ PlmPnoPijk ~ ^PlnPmoPijk ~ ^PloPmnPijk + PklPmoPijn 

+ PklPnoPijm + PkmPloPijn + PkmPnoPijl ~ \pknPloPijm 

— \pknPmoPijl + PkoPlmPijn ~ \pkoPlnPijm ~ \pkoPmnPijl 

2 (pmnoPijkl PijnPklmo ~l~ PjnoPiklm PlnoPijkm PiknPjlmo 
PilnPjkmo PimnPjklo ~t~ Pknopijlrn) ~t~ PinPjklmo PnoPijklm 

These "P-functions have poles of order three for any (n, s)-curve and use seven 
indices. The T, S and V labels were not chosen to signify anything of impor- 
tance. We can of course repeat this procedure to find functions with higher 
orders of poles. Appendix lAl contains details of classes of functions which have 
poles of order four for any (n, s)-curve. Although these functions are not used 
in the remainder of this paper, they have been used in deriving bases for genus 
three curves [15] . 
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Formulation as operators on sigma-functions 



The Hirota operator, introduced in equation 
minant, 

1 1 



, can be expressed as a deter- 



a 

Or, 



d 
diii 



Then one possible generalisation is given by 



T = 



1 


1 


l 


a 


a 


a 


Oil; 


dv , 




a 2 


a 2 


a 2 


0-u 2 


Or 2 


ihi 2 



and the functions 



cr(u) 3 



YiuY jjja(u)a(v)a(w) 



W — V — U 



(29) 



are Abelian with poles of order three and so belong to T(3). Note the similarity 
with the definition of the Q- functions in Definition [5] Although this new class 
of functions is not equivalent to any of the others we introduced above, we do 
have that 



222222 



(tt) = 



1 



<x{u) 



r T 222T 222v{u)a(v)a(w) 



W — V — U 



(30) 



By searching amongst terms with appropriate indices we find that the class can 
be expressed using p-functions as 



1 — ^{^PijPiijj PiijPijj 



2*4) 



= 6Tu 



Hence these S-functions are a subclass of the T-functions introduced above. It 
is anticipated that there may be a more general definition that encompasses 
the entire class of T-functions and indeed that all the functions introduced in 
this section may have alternative definitions involving operators acting on a- 
functions. However, it is not clear yet whether these alternative definitions offer 
any advantage other than brevity of notation. For example, further functions 
in r(3) may be defined by applying further Tm operators to the definition in 
equation (I29[) . However, the resulting functions are a lengthly polynomial in a- 
derivatives involving derivatives of order 12 and it is computationally disadvan- 
tageous to use these in place of the other functions introduced. The derivation 
of Abelian functions using this approach is still a topic under investigation. 



3.3 The trigonal curve of genus four 

In this section we use the generic results of the previous section to derive a basis 
for our example case, the (3,5)-curve. We treat the most general (3,5)-curve, 
namely the curve defined by equation ©. 
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Lemma 5.1 in [3] identified a basis for T(2) in the (3,5)-case as 



{1, pn, pl2, pl3, Pl4, P22, p23, P24, P33, P34, P44, 
<3ll44, Ql244, 02233, Ql444, Q2444}- 



(31) 



Theorem 3. A basis for T(3) associated with the (3,5)-curve is 



{ 



GJJ 



pin, 
P112, 
P113, 

P114, 
P122, 
P123, 
Pl24, 
Pl33, 
P134, 
Pl44, 



P222, 
P223, 
P224, 
P233, 
P234, 
P244, 
P333, 
P334, 
P344, 
P444, 



94Q2444, 
33Q2444, 
92Q2444, 
5lQ2444, 

93<5l444, 
C?2 Q 1444, 
9lQl444, 
94Q2233, 
92Q2233, 



5lQ2233, 
94Ql244, 
(?3<3l244, 
92<3l244, 
<9lQl244, 
94Qll44, 
93Qll44, 
92Qll44, 
5lQll44, 



P 
P 
P 
P 
P 1 
P 1 
P 1 
P 1 
P 1 
P 1 



[3434] 
[2434] 
[2334] 
[2424] 
[2324] 
[1434] 
[1334] 
[2323] 
[1424] 
[1324] 



P 
P 
P 
P 
P 
P 
P 
P 
P 
P 



[1423] 
,[1323] 
,[1414] 
,[1224] 
,[1314] 
,[1223] 
[1313] 
[1214] 
[1213] 
,[1212] 




> . 



Here iSl]) refers to the 16 elements in the basis for T(2) presented above and 
dif refers to the derivative with respect to Ui of the function f . The functions 
p^ kl \Tijkimn,'Pijkimno were all defined previously in this section. 

Proof. The dimension of the space is 3 ff = 3 4 = 81 by the Ricmann-Roch 
theorem for Abelian varieties. All the selected elements belong to the space as 
they are either elements with poles of order two, derivatives of these elements 
or functions explicitly constructed to be in this space. We can easily check their 
linear independence using the cr-function expansion. (Maple is helpful for such 
a computation). 

To actually construct the basis we started by including the 16 functions from 
basis (|3"Tj) for the functions with poles of order at most two. We then know that 
the remaining entries must have poles of order three. We start by looking for 
entries from the set of derivatives of the basis ([31]) . Note that these functions 
need not be linearly independent. In fact in this case we see that four of the 
Q-derivatives have weight —14; 



but that one of these may be expressed using the other three. We are free to 
choose which one of the four we omit from the basis. We test at decreasing 
weight levels and look to see whether these functions can be written as a linear 
combination upon substitution of the series expansions. (We use Maple for these 
calculations and a more detailed discussion of such calculations is given in fTS]. - ) 
Note that while this theorem holds for the general (3,5)-curve in equation 
we need only use the series expansions associated to the cyclic (3,5)-curve in 
equation (| 10[) . This is because if an element cannot be expressed using the basis 
with the restriction on the parameters, then neither will it be expressible with 
the wider set of parameters. Further, we only need to use sufficient expansion to 



S1Q2444, <9 2 Qi444, 93Q2233, ^4<5l244, 
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give non-zero evaluations of the functions considered in order to check whether 
they are linearly independent. 

After examining all these functions we find that 55 basis elements have been 
identified. The next class of functions to be examined were the p[ yfe 'J -functions 
given in Definition [5] We have to test whether they are linearly independent 
both with themselves and with the functions already added to the basis. Another 
simplification that may be made to such calculations is to note that all these 
functions have definite parity according to the number of indices. That is, the 
functions with an odd number of indices are odd and those with an even number 
of indices are even. This can be concluded from the parity properties of the p- 
functions and that all the functions here are constructed from these. So for 
example, when testing to see whether the pW fe '] may be added to the basis we 
need only see whether they are linearly independent of the entries of T(2), since 
all those functions were even and their derivatives odd. We identify a further 
20 basis entries from the p[ Ijfc 'l -functions. Unlike the (3,4)-case, these functions 
were not sufficient to complete the basis. 

To find the final functions we considered some of the new generic functions 
defined in the previous section. First the i3-functions were examined, but in 
this case they were all linearly dependent on the existing basis functions. We 
note that it is possible to reformulate T(3) using 19 of the S-functions in place 
of the 19 Q-derivatives. This formulation may seem less logical but it has some 
advantages as discussed in Section 14.31 To fill the remaining six basis entries 
it is necessary to consider functions with more indices. First the T-functions 
were considered and five more entries identified. We note that the T-functions 
chosen in the basis are not unique. For example, at weight —15 there are 10 
possible distinct T-functions and any one of them may be chosen as the basis 
entry. We choose the function T33344 at random. 

The final entry was then filled by one of the 'P-functions of weight —22. 
Again, this function was just one of a group of 7-index functions, any one of 
which would have been satisfactory. We note that it would have not been 
possible to replace this with a 6-index function however, since such a function 
would be even while this final basis function is odd. This is why none of the 
iS-functions were required in addition to the T-functions. 

□ 

4 Differential Equations 

The Kleinian p-functions satisfy a variety of differential equations which we re- 
view in this section. Some of these differential equations can be found occurring 
naturally in areas of mathematical physics, so the p-functions can be used to 
give solutions to a variety of important problems. In this section we consider 
three main classes of differential equations and present complete explicit sets for 
the functions associated with the (3,5)-curve. We present the relations associ- 
ated with the cyclic restriction of the curve in equation (| 10[) . Relations for the 
general curves can be derived in a similar fashion at a greater computational 
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cost. The sets of differential equations in this paper are presented in decreasing 
weight order as indicated by the bold number in brackets. 



4.1 Four-index relations 

We first consider a set of relations to express the 4-index p-functions. These 4~ 
index relations are the generalisation of equation (|4]) from the elliptic case. We 
aim to express each 4-index p- function as a degree two polynomial in the 2-index 
p- functions, in comparison with equation Then, through differentiation 

and manipulation of this set, we could express all higher index p-functions as 
polynomials in the 2 and 3-index p-functions. Examples of such sets are given 
for the (2,5)-case by Baker in [4] and for the (2,7)-case by Buchstaber, Enolskii 
and Ley kin in [7]. 

There are various ways to derive such relations, with a recent survey given in 
pTj . A constructive way is to consider the basis for T(2) and the set of 4-index 
Q-functions. Each Q-function has poles of order at most two and so belongs 
to the vector space T(2). Hence each Q-function can be expressed as a linear 
combination of basis entries. (The explicit linear combination can be identified 
using the c-expansion as discussed in |18j). 

In the hyperelliptic cases the bases for T(2) can be completed using 2-index 
p-functions and quadratic polynomials constructed of these, (such as the A- 
function discussed in Section 13. So each of the 4-index Q-functions, and 
hence 4-index p-functions, can be expressed as desired giving a set of 4-index 
relations. However, in non-hyperelliptic cases the basis for T(2) can not be given 
purely in 2-index p-functions and so such a set is unobtainable. As discussed in 
the previous section, we can always complete the basis using the Q-functions. 
However, this means that some 4-index Q-functions will be in the basis and 
hence some 4-index p-functions can not be expressed as desired. 

The ability to construct the 4-index relations using only 2-index p-functions 
appears to be a feature unique to the hyperelliptic cases. A more appropriate 
definition for 4-index relations seems to be a set that expresses all the 4-index 
p-functions using a degree two polynomial in the fundamental basis functions. 
(Note that although the polynomials are of degree 2, the Q-functions used will 
only need to appear linearly.) 

Theorem 4. The 4-index relations for the functions associated with the cyclic 
(3,5)-curve are 



( 
( 
( 
( 
( 
( 



4) 
5) 
6) 
7) 
7) 
8) 



P3334 



P2344 



P2444 



P3444 



P3344 



P4444 



-3p 33 + 6pl 4 

3p24 + 6p 34 p44 

-p23 + 2p 34 A 4 + 2p 33 p44 + 4p| 4 
-Q2444 + 6p 33 p34 
Q2444 + 6p24p44 

4pi4 - p22 + 2p24A 4 + 2p23p44 + 4p24p34 



(32) 
(33) 



20 



( — 8) p3333 = 12pl4 - 3p22 + 6p§ 3 



(—25) p\\i2 = 6pnpi2 + (2A + 2A 4 Ai)Qi444 + 6p 33 A 3 Ao 
— 3Q1244A1 

(—28) pun = 6p n - 3p 33 Ai 2 + 8A 4 A <2i 444 + 12p 33 A2A 
— I2Q1244A0 

with the full set given in Appendix\B^ The relations up to weight —19 were 
presented originally in J2J/, but a full set has not been available until now. 

Proof. As discussed above, these follow from the basis for T(2) given in equa- 
tion (|31[) above. The five Q-functions used in the basis appear linearly in the 
equations. The coefficients were determined using the a-expansion, with more 
details on such calculations available in [TJ] for example. 

□ 

4.2 Quadratic relations 

It is natural to next consider a set of differential equations to generalise equa- 
tion (|3]). Such a set should give expressions for the product of two 3-index 
p-functions and so we refer to them as quadratic 3-index relations. The natu- 
ral generalisation would express each product as a degree three polynomial in 
2-index p-functions. Examples of such generalisations have been found in hy- 
perelliptic cases, first for the genus two curve in [3] and more recently for the 
genus 3 case. This is a (2,7)-curve and the corresponding relations were first 
considered in [7] with a complete set recently derived in |15) . 

However, as with the 4-index relations, this generalisation is only possible 
for the hyperelliptic cases. We have explicitly checked that such relations do 
not exist in a variety of non-hyperelliptic cases and so propose the modified 
definition of quadratic 3-index relations to be a set of differential equations that 
expresses all the products of 3-index p-functions using a degree three polynomial 
in the fundamental basis functions. We derive such a set of equations for the 
cyclic (3,5)-case. 

Theorem 5. Each of the 210 products of pairs of 3-index p-functions associated 
with the cyclic (3,5)-curve may be expressed as a degree three polynomial in the 
fundamental basis functions. Some of the relations are given below with the full 
set available online at \16^ . (Note that the polynomial is of degree three but the 
Q-functions need only appear quadratically, or multiplied by a p-function.) 

(-6) p 2 iM = 4p| 4 - 4p 23 - 4p 33 p44 + p^ 4 + 2p 34 A 4 + \\ - 4A 3 

( — 7) P344P444 = 4p 34 p 44 + 2p 2 4p44 ~ P33P34 ~ A 4 p 33 - §Q 2 444 
(-8) P344 = 4p| 4 p44 + 4p M + 4p 24 p34 + pj 3 

(-8) P334P444 = 2p§ 4 p 4 4 + 2p 33 p| 4 - 4p 14 + 2p 22 - 2p 23 p 44 
-p24p34 — 2p§ 3 - P24A4 + 2A4p34p44 
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( — 9) p 3 34p344 = 2(p| 4 + p44p33p34 ~ Pl3 ~ P23P34 + A 4 p| 4 ) + p24p33 
( — 9) P333P444 = 6p44p33p34 ~ 4A 4 p§ 4 + 2p i3 + 7p 23 p34 + 2p 2 4p33 

+A 4 p23 - 2A 4 p 33 p 44 - 2p^ 4 - 2p 34 A| + 6p 34 A 3 + 2A 2 - 2p 44 Q 2444 

( — 9) p244p444 = -2pl3 + P23P34 ~ 2p 2 4p33 ~ A 4 p 2 3 + 2p 34 A 3 - 2A 2 
+4p 24 p 44 + |p 44 Q 244 4 

(—37) pmpii3 = 4:piipi3 — §AiAoQi 44 4 — |A 4 AiQi 444 — 4A0P13Q1244 
+|A 4 Aopi3<5i444 + 4A 4 A 2 AoQi444 + 2Aop 42 + 2AoQ 2444 + 2Aipnpi 2 
-A?pi3p33 - A 2 Aip 33 - 3A 3 AiA p33 + 4A 2 A p 13 p33 + 2AiQi 244 
-6A pnp 22 + 4A pnpi4 + 4A|A p33 - 6A 2 A Qi 244 + 8A 4 Agp 33 

(—39) piupna = 4pf 1 pi2 - |A pnp 2 333 - 2A 3 Aop? 3 - 6A 3 A^p 34 
-2Xj P12P33 + |AiA p 22 33 - T A i A opl 3 ~ 4A^p 23 p 34 + 4A|A p 23 
+4Aop 24 p33 + 8A pi 2 p 44 - A 2 Aip 23 + 2A 3 Aipi3 - 2Aipnpi 244 
-16A 4 AqP 2 3 - 4Aopi 2 pi 244 + 8A3A0P11P33 + 4A 4 AsAopi 2 p33 
+8A 4 A pi 2 p 2 3p33 + 2A 4 A 3 Aipnp33 + 4A 4 Aipnp 23 p 3 3 + X\ P13P23 
+^A 4 A 4 Aopi3 — 6A 3 A 2 Aopi3 + ^A 2 AiA p3 4 + 4A 2 Aopi 2 p33 
-4A 1 A pi 4 p 33 - 5A 3 AiA p 23 + 4AiA pi3p 34 + §AiA p 22 p 3 3 
+4A 1 p 11 p 12 p 44 + 16A pi 2 pi 4 p 24 + 4A pnp 23 p33 - |A 4 Aipiip 23 33 
+8Aipnpi 4 p 24 + ^A 4 A 2 A 4 A - 2Af p 34 + 8Appi 3 - §A 4 A pi 2 p 2 333 
— 18A 4 A3Ag — 2A 4 A^ — ^A^Ao + 8A 2 Aq 

(-42) p\ tl = 4p? x + 8A^p 22 p 33 - 12pnA Qi 2 44 + 8A 2 A§p 34 

-42A 3 Agp 23 - 4A?piip 33 - 4A2A0P13 + 2A 2 A?pi 3 + 8A^pi 3 p 3 4 
-4A?A p34 - 8A^pi 4 p 33 - 8A^A pi3 + 8A 4 A^pi 3 + Af p? 3 - 8A§p| 3 
+16A 2 A pnp33 + 6A3A1A0P13 - 4Afp 23 + I6A2A1A0P23 + 2AgQ 22 33 
-4AiA pi 2 p33 + 4AiA pi3p23 + 8A 4 A pnQi 444 - 4A 4 Af A 
+8A 4 A 2 A2 + A|Af - 4A 3 A? - 4A^A - 27A|Ag + 4AiAg + 18A 3 A 2 A 1 A 

Some of these relations were presented in [5] however the full set and the 
possible structure has not been determined until now. 

Proof. Once again, these relations can be derived through a variety of meth- 
ods as discussed in jlll . They can again be found constructively using the 
er-expansion, although in this case there is no simple linear algebra result to 
dictate that such relations exist, although such a result would surly follow from 
considering a basis for T(6). Instead we may simply search for them using ar- 
bitrary polynomials of the functions from the basis for T(2). The computations 
involved can be heavy and so simplifications are made by ensuring the poly- 
nomials are homogeneous in weight. It is often possible to save computation 
by allocating the cubic terms to cancel higher order poles using Definition 21 
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Another method that minimises computation is to try and find the relations 
through manipulation of the bilinear relations discussed in the next section, (by 
differentiating them or multiplying them by a 3-index p- function). However, it 
is not possible to find all the relations without using the er-expansion. 



In the hyperelliptic cases it is possible to write the set of quadratic relations 
concisely in the form a matrix equation, (see [7] and |15|). Such equations can be 
seen to follow from the representation theoretic approach to Abelian functions 
developed by Athorne in [T] . The corresponding results do not hold for the non- 
hyperelliptic cases but it is possible that similar results may be possible, perhaps 
using a family of matrix equations. This is the subject of current research. 

4.3 Bilinear relations 

The final set of differential equations that we consider here are a set bilinear in 
the 2 and 3-index p-functions. Due to the parity properties of the p-functions 
we know that these bilinear relations cannot contain any constant terms, or 
terms dependent only on the 2-index p-functions. There is no analogue of these 
relations in the genus one case. They have been considered in a variety of cases 



The simplest way to construct these relations is through cross multiplication 
of the 4-index relations. For example, substituting using equations (|32[) and 
([33]) into 



gives the first relation in Theorem [5] below. This was the approach that was 
taken in the papers cited above. However, the existence of Q-functions in the 
4-index relations in the non-hyperelliptic cases means that more care has to be 
taken in the choice of cross products. In higher genus trigonal cases, (or in the 
case where n > 3) the inclusion of further Q-functions in the basis makes this 
method increasingly cumbersome. 

An alternative method to systematically find bilinear relations has been de- 
veloped making use of the B- functions, one of the new generic classes of functions 
discussed in Section [31 These functions are defined in equation (|28p as those 
linear combinations of {pijpimn} which have poles of order at most three for 
a general curve. We can search for bilinear relations as linear combinations of 
£>-functions along with the py^-functions. This approach requires more com- 
putation than taking cross products of 4-index relations. However, as discussed 
in the proof below, it has the advantage of being systematic and allowing us to 
conclude that we have found a complete set of such relations. 

Theorem 6. Every bilinear relation associated with the cyclic (3,5)-curve may 
be given as a linear combination of 65 relations starting with those below. The 



□ 



m 0, Ha, 0, nsi, m and m 
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full set is available online at 

( — 6) = -5P244 - |p333 - P34p444 + P44p344 

( — 7) = 2p234 - A4P344 - P33P444 ~ P34p344 + 2p44p 33 4 

( — 8) = 3A4P334 - 3 p233 - §P33P344 + |p34p334 + |p44p333 

— P24P444 + P44P244 
( — 9) = -2pi44 - P33P334 + P34p333 + p24p344 ~ P34p244 
( — 9) = 2p224 - 2pi44 - A 4 p244 - P23P444 + P34p244 

+ 2p44p234 - 2p24p344 

( — 32) = pl2plll - PllPll2 + -|A4Ao(2A4Aopi3p334 - 2p34pi 33 
-2pl2p444 - 2p24pl44 + ^P±Ap\2A) + 5A3A1 (pl3p334 - P34P133) 
+ A|A (|p44p224 - |p22p444 ~ |p24p244 ~ P23P334 + P34P233) 
+ gA^A (2A|A p 3 4p334 + 2A|A p44p333 - 4p33p344 ~ A|A p233) 
+ 3 1 Ao(4pl4pl44 - 4p44pi22 — 4p 2 4pl24 + 2p 2 2pl44 + 5p 2 3pl33 
+5pl2p333 - 8p44pll4 - 10p33pl23 + 6pi2p244 + 4pnp44 4 ) 
+A3A (p22p444 + P24P244 ~ 2p4 4 p 2 24 + 3p 2 3p334 — 3p 3 4p 2 33 
+ 5P133) + |A?(2p 233 - P34P334 + 2p 33 p3 4 4 ~ P44P333) 
+ iA 4 A 3 A (p233 - 2p 3 4p334 ~ 2p44p3 33 + 4p33p 3 4 4 ) 
+ ^Ai(pi2p224 — P24P122) — |A|A3Aop334 + gA4A^p334 
+2AgA p334 + |A|A p334 + A1A0P334 — 2A 2 A p233 — 5A2A1P133 

A number of these relations were originally derived in [5] but this complete set 
was not available until now. 

Proof. These relations were derived alongside an alternative basis for r(3), con- 
structed as follows: We start by including the 3-index p-functions in the basis. 
We then consider decreasing weights in turn and look to see how many B- 
functions at that weight may be included in the basis. Note that unlike the 
p-functions the order of the indices in the B- function is relevant. The definition 
of the £>-functions can simplify for specific indices and so we first test to see 
which £>-functions at a given weight are distinct. We then test their linear de- 
pendence with themselves and with the pijk ■ Note that when doing this we need 
not consider the entries in T(2). Although they are included in the basis they 
are all even functions and so cannot be linearly dependent on the S-functions 
which are all odd. We test the linear independence using the cr-expansion as 
discussed in the proof of Theorem [3] It is important to make use of weight 
simplifications when expanding products of series to reduce computation. 

We find that the following 19 functions may be added to the basis, at the 
weights indicated. Note that they play the role of the 19 Q-derivatives in the 
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alternative basis in Theorem [3] 



( — 8) $33344 
(-9) #23444 



(-16) B12334 

( — 17) #12333, #12244, #11444 



— 11) #23334, #22444 



( — 18) -B11344 

(-19) £12233 



— 12) #13444 

— 13) #13344 




— 14) #22244, #13334, #12444 

— 15) #12344 



Now, any bilinear relation between the 2 and 3-index p-functions must have 
the order 4 and 5 poles in the pij pki m -terms canceling. These terms must be 
hence composed of £>-functions: The only alternative would be if there were 
combinations with reduced orders of poles specific to the curve. It is simple to 
check that no such combinations exist in the (3,5)-case as we can check using 
just the leading order terms of the cr-expansion, given in (|22p . So all bilinear 
relations are given by linear combinations of B- functions and pijk- There will 
be one such relation for each of the distinct B-functions not in the basis, and 
every other bilinear relation will be a combination of these. The actual relations 
were derived as a by-product of the derivation of the alternative basis to T(3). 



In general, this approach allows us to derive all independent bilinear relations 
and conclude we have the full set. However, care must be taken in drawing 
the conclusion that all bilinear relations have pole cancelation in the form of 
£>-functions. As discussed in the proof, we need to check that their are no 
cancelations specific to the curve. In the (2,5)-case for example, there exists 
the function A in equation (|27l) which has poles of order 3 in general but poles 
of order 2 on the (2,5)-curve. Hence the derivatives of A are sums of pijpkim 
which are not £>-functions and yet have poles of order three. This leads to extra 
bilinear relations. The existence of functions like A and the following extra 
bilinear relations appears to be a feature unique to the hyperelliptic cases. 

We note that there are no bilinear relations at weights —5 or —35 in associ- 
ation with Z?44444 and #nm as may be expected. This is because the definition 
of the S-functions collapses to zero in the case where all the indices are the 
same. This also explains why there can be no bilinear relations in the elliptic 
(genus one) case. 

5 Addition Formulae 

Here we discuss the addition formulae satisfied by the Abelian functions and 
present a new formula associated with the cyclic trigonal curve of genus four. 
We start by considering the formulae which generalise equation ([5]) from the 
elliptic case. Such a formula will express the ratio of functions 



□ 



a(u + v)a(u — v) 



a(u) 2 a(v) 2 
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as a sum of products of pairs of functions from the basis for T(2), one a function 
of u and the other of v. This must be so since the ratio above can be shown to 
be Abelian in both u and v and it has the correct pole structures, (see [TS] for 
details). We can go further and use the parity property of a(u) to check that 
the function is either symmetric or anti-symmetric with respect to the change 
of variables (u,v) i-> (v,u) when the tr-function is odd or even respectively. 

The coefficients in the polynomial can be explicitly determined using the a- 
expansion. (See [TH] for details of such calculations). Such formulae have been 
found for a variety of cases in g], [7J, [5], [IB], J7J. In particular, the formula 
for the cyclic (3,5)-curve was derived in Section 8 of [5] and is given by 

a(u + v)a(u - v) 

= f(u, v) + f(v, u) 

where 

f(u, v) = p4i(u)p n (v) + p22(u)p u {v) - pi 2 (u)p 2i (v) 

-2p 14 (u)p 14 (v) - \Ql22i{u) + Ql244(v)p34(u) 

+ ^Q2233(u)p33(v) - ^Q2444(v)pl 3 (u) - \ Q 2 333 («) P23 (v ) 

+ lQ2333{v)p34{u)\4 - |P13(m)P33(«)A 4 ~ |P33(m)P34(^)A 2 

-lQ2333(w)A 3 ~ P33(m)P34('u)A4A3 + ( |A X - 5A4A2 + |A§^P 33 (m). 

Note that the (3,5) c-function is even and hence the addition formula here 
is symmetric in (u,v). (The reason for the difference between this formula 
and the one is [5] is that some formulae in |S] were calculated using different 
Q-functions in the basis. Our paper makes choices consistent with the basis 
derived in Section 5 of [5].) 

In cyclic non-hyperelliptic cases there are more addition formulae associated 
with the functions, resulting from automorphisms of the curve equation. Such 
addition formulae were the topic of |14j which gave a thorough treatment of the 
genus one and two cases. We will present a new genus four addition formulae 
associated with the cyclic (3,5)-curve. It is related to the automorphism of the 
curve (jlOp given by the operator 

/ 2ni 

[C] : (x, y) i-> (x, (y), where C = exp f — 

So C is a cube root of unity and [£] an operator which multiplies y by the root 
leaving the curve unchanged. We extend this notation to define the sequence of 
operators and automorphisms, 

[?}:(x,y)^(x,(iy), for j 6 Z. 



We can check using the basis of differentials (1121 that these operators act on 
the variables u as follows. 

[C> = {( j ui, Cu 2 , C 2l u 3 , ( J u 4 ). (34) 
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The action of such operators on the lattice A is stable, (it moves the points 
around but does not change the overall lattice) . This can be checked by consid- 
ering the effect on the individual elements of the period matrices, (see [21] for 
more details). We can now derive the following result for the a- function which 
follows Lemma 4.2.5 in [2~i] . 

Lemma 7. The a-function associated to the cyclic (3,5)-curve satisfies 

a([C j ]u) = C 2j cT(u). (35) 

Proof. Consider the quasi-periodicity of a"([£ J ']ti). If £ is a point on the lattice 
then 

a([C j ](u + £))=a([C j ]u+[^]£). 

Since the lattice is stable under the action we know that [C]£ is also on the 
lattice. Hence by (fn 



a([C}(u + lj) = X (K j }t)o-([(i]u)e W (u + -) , [C]/)~ 

In |24j the author shows that for an automorphism of a cyclic curve we have 

L([Qu,v) = L(u, [C» 

and hence 

L([C>, [C» = L(u, iC j ][C j }v) = L(u,v). 
Therefore, we have 

a(lC}(u + £)) = x([CW[C>) exp 
We now consider the quotient 



L(u+-,£ 



a((u + £)) xW <u) a(u) 

since = ±1. So we see that the function 

^ac>) 

o-{u) 

is bounded and entire (since the zero sets coincide). Hence, by Liouville's the- 
orem, the function is a constant. Using the Schur-Weierstrass polynomial (|22|) 
we see that this constant is C, 2 ^ . 

□ 

We can now derive the addition formula associated with these automor- 
phisms. 
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Theorem 8. The functions associated to the cyclic (3,5)-curve satisfy 



where 

f( U , V) = [p 24 + P 2 1 + Pl8 + Pl5 + Pl2 +P 9 +P 6 +P 3 + P ] («, V) 

and the polynomials Pk(u,v) are as presented in Avvendix[C\ 

Proof. Denote the left hand side of the formula by LHS(m,i>). Using Lemma 
[7] and the parity property of the c-function we first check that LHS(u,v) is 
symmetric under (it, v) h-> (v, u). Next consider the affect of u (-» u + £. 

IHW v o-(u + e + v)a(u + i + [t]v)a{u + £ + [( 2 }v) 
Ltib(u +l,v) = 



LHS(«,«) 



a(u + £) 3 a(v) 3 
x (e)e^^ u + v + e /^x(£)e^ [u+[c]v+e/2] x{i)e^ [u+[C ' 2]v+e/2] 



x (£)3 e 3*[«+f/2] 

= LRS(u,v)e^ 1+ ^ + ^. 
However, from (l34l) we see that 



Vl (i + c + C 2 )\ /o 

^(i + [c] + ic 2 ])= I v 2 (i+c+e) = 
«3(i+c 2 + o; \o 

and so LHS(u, v) is Abelian with respect to u. Further, since it is symmetric 
in (u,v) we can conclude that it is Abelian in v also. It is now clear that 
LHS(m, v) may be expressed as 

LHS(u,t;) =^c. i A l (M)B l («) 

i 

where the Cj are constants and the functions Ai(u),Bi(v) belong to the basis 
for F(3), presented earlier in Theorem [3] 

To determine the constants c, we use the a-expansion. The computations 
involved can be heavy and so it is essential that we take into account all the 
available simplifications. We have already noted that LHS(w, v) is symmetric in 
(u, v), but we can also check that it is even with respect to (u, v) m> (— u, —v). 
Hence we need only consider combinations of functions that preserve these two 
properties. Further, we can check that it has total weight —24 and so only 
consider those combinations with this weight (including when combined with 
monomials in the curve constants). Together these simplifications drastically 
reduces the number of possible terms. To further ease the time and memory 
constraints we implement code in Maple to efficiently expand the products of 
series so that only the relevant terms are considered. 
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We find that LHS(m,i>) is given as stated. For simplicity we group together 
the terms with common weight ratios into the polynomials Pk{u, v) which con- 
tain the terms with weight — k in the Abelian functions and weight — (24 + k) 
in the curve parameters. The polynomials are presented in Appendix [C] 

□ 

This formula is actually a reduction of a more general formula in three sets of 
variables expressing 

a(u + v + w)a{u + [Qv + [C 2 ]w)a{u + [( 2 }v + [Qw) 
a(u) 3 a(v) 3 a(w) 3 

as a sum of products of triplets of functions from r(3), one each in u,v and 
w. The existence of these formulae were predicted in [5] but we can only now 
derive them since they require the basis for r(3), calculated in Section [3] using 
the new classes of Abelian functions that are the subject of this paper. 

We finish by noting that there are further addition formula associated with 
further reductions of the (3,5)-curve. The curve y 3 = x 5 + Ao will have a family 
of automorphisms given by 

[i J ] : (x,y) i ^ (rfx,y) 

where 77 is a fifth root of unity. We can check that the ratio 

a{u + v)a(u + [rj]v)a(u + [f] 2 ]v)a(u + [rf]v)a{u + [i] 4 ]v) 
<j(u) 5 (j(v) 5 

is Abelian in both u and v and so a similar approach could be employed to 
derive the formula, following the derivation of a basis for T(5). However both 
these computations would be computationally intensive. 

A Generic Abelian functions with poles of order 
four 

In this Appendix we repeat the procedure described in Section [3T21 to calculate 
Abelian function which have poles of order four for any (n, s)-curve. 

The first non trivial case occurs when we allow five indices. We may con- 
struct polynomials from 5-index p-functions and the products of a 2 and 3-index 
p-function. There is only one independent non-zero combination which has poles 
of order four, given by 

Pijklm + PijPklm ~ ^SplmPijk- 

(We note that the coefficient 13 may seem unexpected but it can be easily 
verified that it is necessary to balance the poles of order 5 in the other two 
terms.) 
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With six indices we can build polynomials from the set of functions 

{pijklmm PijkPlmn: Pij Pklmn 7 PijPklPmn\ • 

There are three independent non-zero combinations which have poles of order 
4. These are given by 

Pijklmn "t - pij pklmn 2 (^(pjkpilmn "t - Pjlpikmn "t - PjmPikln 

+ PjnPiklm) + 5(pfcZ Pijmn + Pkmpijln ~\~ PknPijlm 4" Plmpijkn 
+ PlnPijkm + PmnPijkl)), 

Pij Pklmn j(3pijkplmn + 3pijl 

Pkmn H~ ^Pijmpkln H~ ^Pijnpklm 
PiklPjmri PikmPjhi PiknPjlm PilmPjkn Pilnpjkm PimnPjkl) i 
PijPklpmn ~ \{pijkplmn + Pijlpk mn ~\~ Pijmpkln H~ Pijnpklm Piklpjmn 
PikmPjln PiknPjlm PilmPjkn PilnPjkrn PimnPjkl) • 

Note that we take linear combinations of these to define the class, 

Fijklmn = PijPklp mn PijPknPlm PilPjkPmn "r PilPjnPkm 

+ PimPjlPkn — PimPjnPkl + PinPjkPlm ~ PinPjlPkm, 

which were used to complete the basis for T(4) in the (3,4)-case as detailed in 
|15j . This procedure may be continued indefinitely. We have compiled functions 
with up to 8 indices and these can be found online at [16 . One class of particular 
note is 

Gijklmnop — Pijklmnop ^(pijnopklmp "t - PijnpPklmo ~t~ PijlpPkmno 

+ Pij mnpklop H~ Pijmopklnp ~t~ Pijmppklno Pijkmplnop T" Pijknplmop 
~^~Pijkoplmnp PijkpPlmno H~ Pijlmpknop H~~ Pijlnpkmop H~ Pijlopkmnp 
\~PijklPmnop H~~ PijopPklmn H~ Piklmpjnop ~\~ PiklnPjmop H~ Piklopjmnp 
+ PiklpPj mno H~ Pikmnpjlop H~ Pikmopjlnp H~~ PikmpPjlno T" piknopjlmp 
+ PiknpPjl mo ~T PikopPjlmn PilmnPjkop H~~ PilrnoPjknp H~ PilmpPjkno 
+p ilnopjkmp H~~ PilnpPjkmo \ PilopPjkmn H~~ PimnoPjklp H - PimnpPjklo 
~^~PimopPjkl7i PinopPjklm) ■ 

which was nsed to complete the basis for T(4) in the (2,7)-case as detailed in 
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B The 4- index relations 



The full set of 4-index relations for the functions associated with the cyclic 
(3,5)-curve are given below. 

( — 4) P4444 = 6p| 4 - 3p 33 

( — 5) P3444 = 6p 3 4p44 + 3p 2 4 

( — 6) P3344 = 4p34 - p23 + 2p 34 A 4 + 2p 33 p 4 4 

( — 7) P3334 = 6p33p 34 - ^2444 

( — 7) P2444 = 6p24p44 + Q2444 

( — 8) p2344 = 4pi4 - p22 + 2p 24 A 4 + 2p2 3 p44 + 4p24p34 

(-8) p3333 = 12pl4 - 3p 22 + 6p§ 3 

( — 9) P2334 = 2pi 3 + 3p 34 A 3 + A 2 + 4p 23 p34 + 2p 24 p33 

( — 10) p2244 = IQ2444A4 + P33A3 + 2p22p44 + 4p24 + §Ql444 

( — 10) p 1444 = Q1444 + 6pi4p 44 

( — 10) P2333 = 6p 23 p33 - 2Q 144 4 + 3p 33 A 3 

( — 11) p2234 = 4p23p24 ~ 2pi2 + 4p M A 4 + 3p24A 3 - 2p 44 A 2 + 2p22p34 
( — 11) P1344 = 4pi4p 34 - pi2 + 2pi 4 A 4 + 2pi3p 44 

( — 12) P1334 = 2pi 3 A 4 + f P23A3 + 2p 34 A 2 + A 4 A 2 + 4p 13 p 34 

+ 2pl4p33 - 5Q2233 
( — 12) p2233 = Q2233 + 2p 22 p33 + 4p| 3 
( — 13) pl333 = 6pi3p33 - 2A4Q1444 + 3Qi244 
( — 13) pl244 = Ql244 + 2pi2p44 + 4pi4p24 

( — 13) p2224 = 4A 4 (3i444 + Q2444A3 + 6p 33 A 2 - 6Q1244 + 6p22p24 

( — 14) P2223 = 6P22P23 - 6pn + 6pi 4 A 3 + 3p22A 3 - 6p 44 Ai 

( — 14) pi234 = 2pi4p23 - pll + 3pi4A 3 ~ P44A1 + 2pi2p34 + 2pi 3 p 2 4 

( — 15) P1233 = 3pi 3 A 3 + 2p 34 Ai + A4A1 - 3A + 2pi 2 p 33 + 4pi 3 p 23 

( — 16) p2222 = 4Q2444A2 — I2Q1244A4 + 8A4 2 Qi444 + 2A3Q1444 

-12p3 3 A 3 2 + 12p33A 4 A 2 + 6p 33 Ai + 6Q1144 + 3Q2333A3 + 6p| 2 
( — 16) P1144 = Q1144 + 2piip 44 + 4p 2 4 

( — 16) pi224 = —Qllii — 5<32333^3 + §P33^3 2 + 3p33Ai 
+4pl2p24 + 2pi4p22 

( — 17) P1223 = -2pnA 4 + 3pi 2 A 3 + 4pi 4 A 2 - 2p 2 4A x - 6p 44 A 

+4pl2p23 + 2pl3p22 

(—17) P1134 = 2pi 4 A 2 - P24A1 + 2piip 34 + 4pi 3 pi4 
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— 18) pn33 = 2pnp33 + 4pi 3 + 2p 13 A 2 - P23A1 + 6p 34 A 

— 19) pH24 = §A 2 Ql444 - 3Q2444A1 + 4p 33 A + 2pnp24 + 4p i2 pl4 

— 19) pi222 = — 3Qi244-^3 + 2A4A3Q1444 + 2Q 2 444^1 + 6P33A4A1 

+ 6p33Ao + 6pi2p22 

—20) P1123 = 2pi 2 A 2 - 6P24A0 - p2 2 Ai + 4pi 4 Ai + 2pnp23 + 4p i2 pi3 

— 22) pn22 = p33^3^1 + 8P33A4A0 + |A 4 A 2 (5 1 444 — 2Q 12 44A 



2 



-2Q2444A0 + 2pnp 22 + 4pf 2 + |AiQ 



^1444 



—22) puu = -2Q2444A0 + 6pnpi4 + A1Q1444 

—23) pni3 = 6pnpi3 - 6P22A0 + 6pi 4 A + 3pi 2 Ai 

—25) pni2 = 2A0Q1444 + 6P33A3A0 + 2A4A1Q1444 + 6pnpi2 - 3Q1244A1 

—28) pun = §p\ x - 3p 33 Ai 2 + 8A4A0Q1444 + 12p3 3 A 2 A - \2Qi 2iA X 

These relations are also available online at [T5] . 

C The addition formula 

The addition formula presented in Theorem [5] is constructed using the following 
polynomials. 

P 24 (u,v) = - ±7222222 («) + §71 14444 (v ) p23 («) + §7l33344 («) Pl3 (u) 
+ |722 2 224(m)P34(v) + g p!22 («) P224 («) + | Pll4 («) P224 («) 
-gPl24(«)pl24(«) - iPl44(«)plu(u) - g p 3 33 («)pll2 («) 
+ iPl22(M)pl44('w) - |P222(W)P124('W) ~ | Pll2 («) P244 0) 
+ Jp444(u)pm{v) + JgQ 2233 {u)Q 2233 (v) ~ ^d 3 Q 124 4 (u)8 3 Q 2 444 (v) 

-p^ 13 Hv)p 23 (u) + \p^\v)p^\u) + \p^\v)p 13 {u) 

+ l p [3434] (t> j p [1313]( u ) _ l p [2323] (M)p [1334] (t;) _ I p [1334] („) p [1334] („) 

-|^3<9i144(^)P333(m) - ^pl22 (m)9 3 Q2444 («) + 2 ^P2 2 2(m)9 3 Qi444('w) 
+ !<9 3 Qi244(»P 2 24(m) ~ ^d 3 Q 12 4 4 (u)p 14i (v) + ^3 Q 2 444 («) Pll4 ) 
-f2"93Qi444(M)9 3 (5l444('y) ~ gPl24(«) ,9 3Ql444('") 

P al («,t;) = (ip 13 (M)p[ 1334 ](«) ~ ipl3(«)Q223s(«) -P34NP [13131 W 
- 1 ^P2 2 4(m)9 3 Qi444('w) + \p!44{v)d 3 Qi444{u) + \ pn 4 («) p 333 (v ) 
+ ^9 3 Q 24 44('w)93(5i444(m) + §7222224 («) ) A 4 
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Pia{u,v) = ( |7Il4444(») - |7l33344(v)P34(m) - §P333(m)P124(« ) 
-|p 23 («)Q2233(«) ~ P^(V) + |p[1323] (M)p34(v) 

+ §p 13 (u)pi 3 (u) + 3jPa33(w)9jQi444(u))A3 - ip[ 1313 ]( M )Ai 

Pl6(u,w) - (|p [1323] («) + |P333(«)P144(m) + §Q 2233 (v)P3 4 (m) 

-P34(M)P [1334] («) " |P34(«)P [23231 (^) + |p23(«)pi3(t;) 
+ ip[ 3434 ](«)pi 3 (u) + |ri33344(«))A 2 

+ ( |p 23 (ti)pl3(t») - fTl33344(«) - |P [13231 N 
-|P34(«)P [13341 (^) " |p34(tt)Q 2a 33(t>))A 4 A3 

P 12 (u,v) = ±Q 2233 (u)A 2 + \p^ 2 \v)\ x - iQ 2233 («)A 1 

+ |P34(m)P [23341 ('w)Ai - \p-2^(u)p 333 (v)\ x + f p 3 4(w)pl3 

+^p 23 (M)p 23 («)A 2 - ip 23 (M)p[ 3434 ](^)A 1 + iQ 22 33(«)A 4 A 2 
-lp[ 1334 l(«)A|A 3 - ip[ 2323 l(«)A 4 A 2 - |p[ 2334 ](^)p 34 ( M )A 2 
- 3 <W«)A|A 3 + ip[ 1334 l(M)A 4 A 2 - |p 13 («)p 34 («)A 2 
-\p^\u)\\ - ip[ 2323 l(«)A| - ip 3 4^)pi3(«)A 4 A 2 

P 9 (U,V) = (jp3±(v)p23{u) + jp 333 (v)p 4ii (u) - |pl 3 (w))A 

( ipl 2334 l (u) + 3p 13 (v) - fpl 3434 ! (v)p 34 (u) - lp 23 (v)p 34 (u) ) A 3 A 2 
(2p ls («) + |p [2334] (u) + fp3 4 (v)p23(M))A 4 Ai - ipi 3 (M)AiA 2 
( \p2 3 {v)p 3i {u) - 2p 13 (u) - fpl 2334 ! (u) ) A 4 A| + p 13 (M)A 3 A 3 

P 6 (u, v) = ( A3A1 - ±A 2 + |A 4 A )p 23 («) + ( f AlA 2 - f A 3 + A 2 A 4 )p 23 (w) 
+p[ 3434 l(«)A 2 - ip[ 3434 l(«)A 3 A 1 - |p[ 3434 ](«)A 4 A 3 A 2 
+ (|Al - |A|Ai - 3A 4 A + 4A4A3A2 )p 3i {u)p 3i {v) 

P 3 (u,v) = (|A 3 A - |A 2 Ao - ±A§A 2 - §A 3 Ai - ±A 4 A 2 )p 34 (u) 
+ ( f A|A 3 A 2 + |A 4 A 3 Ai )p 34 (u) 

Po(u,v) = — ^A^A 2 , + 3A3A1 + ^ A4 Ao — |A 2 Ao — 1 jA 4 A|A 2 + tjA 3 A3A 2 
— i;A 2 A3Ai + |A 4 A 2 Ai + |A 4 AsAo — 2A 2 

These polynomials are also available online at [IB] . 



33 



References 



[1] C. Athorne. Identities for hyperelliptic p-functions of genus one, two and 
three in covariant form. Journal of Physics A: Mathematical & Theoretical, 
41:415202, 2008. 

[2] H. F. Baker. Abelian functions: Abel's theorem and the allied theory of 
theta functions. Cambridge University Press, 1897 (Reprinted in 1995). 

[3] H. F. Baker. On a system of differential equations leading to periodic 
functions. Acta Mathematica, 27:135-156, 1903. 

[4] H. F. Baker. Multiply periodic functions. Cambridge University Press, 
Cambridge, 1907 (Reprinted in 2007 by Merchant Books. ISBN 193399880). 

[5] S. Baldwin, J. C. Eilbeck, J. Gibbons, and Y. Onishi. Abelian functions 
for cyclic trigonal curves of genus four. Journal of Geometry and Physics, 
58:450-467, 2008. 

[6] S. Baldwin and J. Gibbons. Genus 4 trigonal reduction of the Benny equa- 
tions. Journal of Physics A: Mathematical and Theoretical, 39:3607-3639, 
2006. 

[7] V. M. Buchstaber, V. Z. Enolskii, and D. V. Leykin. Klcinian functions, 
hyperelliptic Jacobians and applications. Reviews in Math, and Math. 
Physics, 10:1-125, 1997. 

[8] V. M. Buchstaber, V. Z. Enolskii, and D. V. Leykin. Rational analogs 
of Abelian functions. Functional Analysis and its Applications, 33:83-94, 
1999. 

[9] V. M. Buchstaber, V. Z. Enolskii, and D. V. Leykin. Uniformization of 
Jacobi varieties of trigonal curves and nonlinear equations. Functional 
Analysis and its Applications, 34:159-171, 2000. 

[10] K. Cho and A. Nakayashiki. Differential structure of Abelian functions. 
International Journal of Mathematics, 19:145-171, 2008. 

[11] J. C. Eilbeck, V. Z. Enolski, and J. Gibbons. Sigma, tau and abelian 
functions of algebraic curves. J. Phys. A: Math. Theor., 43:455216, 2010. 

[12] J. C. Eilbeck, V. Z. Enolski, S. Matsutani, Y. Onishi, and E. Previato. 
Abelian functions for trigonal curves of genus three. International Mathe- 
matics Research Notices, page Art. ID: rnml40 (38 pages), 2007. 

[13] J. C. Eilbeck, V. Z. Enolskii, and D. V. Leykin. On the Klcinian con- 
struction of Abelian functions of canonical algebraic curves. In D Levi and 
O Ragnisco, editors, Proceedings of the 1998 SIDE III Conference, 1998: 
Symmetries of Integrable Differences Equations, volume CRMP/25 of CRM 
Proceedings and Lecture Notes, pages 121-138, 2000. 



34 



[14] J. C. Eilbcck, S. Matsutani, and Y. Onishi. Addition formulae for abelian 
functions associated with specialized curves. Phil. Trans. R. Soc. A, 
369(1939):1245-12638, 2011. 

[15] J.C. Eilbeck, M. England, and Y. Onishi. Abelian functions associated 
with genus three algebraic curves, preprint arXiv:1008.0289, 2010. 

[16] M. England. 

http: //www. maths .gla. ac .uk/~mengland/Papers/2011_DBAF/. 

[17] M. England. Higher genus Abelian functions associated with cyclic trigonal 
curves. SIGMA, 6:025, (22pp), 2010. 

[18] M. England and J. C. Eilbeck. Abelian functions associated with a cyclic 
tetragonal curve of genus six. Journal of Physics A: Mathematical and 
Theoretical, 42:095210 (27pp), 2009. 

[19] M. England and J. Gibbons. A genus six cyclic tetragonal reduction of the 
Benney equations. Journal of Physics A: Mathematical and Theoretical, 
42:375202 (27pp), 2009. 

[20] S. Lang. Introduction to algebraic functions and Abelian functions. Num- 
ber 89 in Graduate Texts in Mathematics. Springer- Verlag, 2nd. edition, 
1982. 

[21] A. Nakayashiki. On algebraic expressions of sigma functions for (n, s)- 
curves. Asian Journal of Mathematics, 14(2):175-212, 2010. 

[22] A. Nakayashiki. On hyperelliptic abelian functions of genus 3. Journal of 
Geometry and Physics, 61(6):961-985, 2011. 

[23] A. Nakayashiki and F. Smirnov. Cohomologies of affine hyperelliptic Jacobi 
varieties and integrable systems. Comm. Math. Phys., 217:623-652, 2001. 

[24] Y. Onishi. Complex multiplication formulae for hyperelliptic curves of 
genus three. Tokoyo Journal Mathematics, 21:381-431, 1998. A list of 
corrections is availabe from http://web.cc.iwate-u.ac.jp/~onishi/. 



35 



